An impurity in a Fermi sea on a narrow Feshbach resonance: A variational study of 

the polaronic and dimeronic branches 
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We study the problem of a single impurity of mass M immersed in a Fermi sea of particles of mass 
m. The impurity and the fermions interact through a s-wave narrow Feshbach resonance, so that 
the Feshbach length R, naturally appears in the system. We use simple variational ansatz, limited 
to at most one pair of particle-hole excitations of the Fermi sea and we determine for the polaronic 
and dimeronic branches the phase diagram between absolute ground state, local minimum, thermo- 
dynamically unstable regions (with negative effective mass), and regions of complex energies (with 
negative imaginary part). We also determine the closed channel population which is experimentally 
accessible. Finally we identify a non-trivial weakly attractive limit where analytical results can be 
obtained, in particular for the crossing point between the polaronic and dimeronic energy branches. 

PACS numbers: 03.75.Ss, 34.50.Cx, 05.30.Fk, 71.10.Ca 
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I. INTRODUCTION 

The physics of strongly interacting Fermi gases is 
now experiencing a rapid development. Experimental- 
ists have succeeded in obtaining highly degenerate and 
strongly interacting samples, providing new perspectives 
in the study of Fermi superfluidity, including BCS-BEC 
crossover: Thanks to magnetic Feshbach resonance tech- 
niques one can indeed control at will the interaction 
strength and even reach the unitary limit, where interac- 
tion effects are maximal [l|, . Furthermore very precise 
measurements of the many body properties of these sys- 
tems are at reach: Recently the equation of state of the 
gas at the unitary limit or within the whole BEC-BCS 
crossover was measured with an uncertainty at a few per- 
cent level 

Of particular interest is the case of spin imbalanced 
Fermi gases, that is with a larger number of spin up par- 
ticles (iV-|- ^ iV^). In this spin polarized regime, the usual 
Cooper pairing scenario no longer applies, and this raises 
the question of the possible phases of the system. Apart 
from the FFLO phase 0], and more exotic phases Q, it 
was proposed recently that, the minority atoms dressed 
by the Fermi sea of the majority atoms form a weakly 
interacting gas of quasi-particles called polarons 's', ^] . 

Whereas the basic properties of these polarons, such as 
their binding energy with the Fermi sea, and their effec- 
tive mass, are well understood in the limit where the ma- 
jority atoms interact with the minority atoms resonantly 
on a s-wave wide Feshbach resonance 0-[l3> their prop- 
erties have not been studied yet when the Feshbach res- 
onance becomes narrow. A narrow Feshbach resonance 
is characterized by a length > 0, called the Fesh- 
bach length, which vanishes in the broad resonance limit 
R* = JJi]- Therefore the presence of a positive i?* in- 
troduces a new length in the system, in addition to the 
s-wave scattering length a, which may be of physical im- 
portance depending on how it compares with the other 
physical parameters. It is then natural to ask whether 
the impurity problem on a narrow Feshbach resonance 



presents new physics. Of particular interest is the case 
of fermionic mixtures of ^Li and ^°K, which in the last 
decades of experiments have been extensively realized. It 
turns out that all the Feshbach resonances in such mix- 
tures are narrow [ij, characterized by a value of the 
Feshbach length i?* > lOOnm, which is much larger than 
the typical van der Waals length b ~ 2nm. It is there- 
fore expected that a fermionic mixture of ^Li and "^^K 
will present some physical effects due to the narrowness 
of the Feshbach resonances, as recently demonstrated in 
the experiment of (Tsj . 

In this paper we study the single impurity problem on 
a narrow Feshbach resonance, that is a single impurity 
interacting through a s-wave narrow Feshbach resonance 
with a Fermi sea. As shown in 10] on a broad Fesh- 
bach resonance this problem has two energy branches, 
one that continuously connects to the free impurity in 
the weakly attractive limit (a — > 0^), that we call the 
polaronic branch, and the other one that continuously 
connects to the free impurity-atom dimer for a — > 0"*", 
which we call the dimeronic branch. We find that the 
scattering length value at the crossing point between the 
two energy branches strongly depends on the Feshbach 
length i?*, as well as on the mass ratio between the im- 
purity and one of the fermions. In particular, for all 
the mass ratios considered in this work we find that as 
Rf: increases, the domain of the polaronic ground state 
shrinks and shifts progressively towards the weakly at- 
tractive region. For any negative scattering length, we 
find that for a sufficiently large i?* the ground state is 
dimeronic. Furthermore, we identify a non-trivial weakly 
attractive limit, a — > 0^ with i?*a fixed, which allows us 
to solve in an asymptotic expansion the single impurity 
problem and to obtain the crossing point between the 
polaronic and dimeronic branches. 

The paper is organized as follows. After introducing 
the narrow Feshbach resonance on a two-channel model 
in section |TT1 we discuss the variational ansatz in sec- 
tion mil We present the physical motivation behind the 
ansatz, we show that the problem falls naturally into the 



2 



well studied class of a discrete state coupled to a con- 
tinuum, and we derive the variational eigenvalue equa- 
tions. In section IIVI we numerically solve the eigen- 
value equations and we calculate the properties of the 
system, including the crossing point, the effective mass 
and the closed-channel molecule population. In sectionlVl 
we develop a non-trivial weakly attractive limit that we 
can treat analytically in a perturbative way. We then 
show that the analytical solution compare remarkably 
well with the numerical results of section IIV[ even for 
parameter values that are not in the asymptotic region. 
We conclude in section IVll 



II. TWO-CHANNEL MODEL HAMILTONIAN 

We consider a Fermi gas of same-spin-state particles of 
mass m perturbed by the presence of an impurity, that is 
an extra, distinguishable particle of mass M. We assume 
that there is no interaction among the fermions, and that 
the impurity interacts with each fermion resonantly on 
a s-wave narrow Feshbach resonance. We use an usual 
two-channel model to describe this resonant interaction 
[l6l - l25j : The particles exist either in the form of atoms in 
the so-called open channel (annihilation operators u and 
d for the fermions and the impurity, respectively), or in 
the form of a tightly bound molecule of a fermion with the 
impurity (annihilation operator b) , in the closed channel. 
We shall take the limit of an infinitely narrow Feshbach 
resonance jl^, neglecting the direct interaction between 
the atoms in the open channel, and including only the co- 
herent interconversion of the closed-channel molecule into 
a pair of atoms. The resulting second-quantized Hamil- 
tonian for the system enclosed in a cubic quantization 
volume V with periodic boundary conditions is then 
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-== E^(^^<=l)(^k+k'"k4' +h.C.) , (1) 
k,k' 



where wj^. creates a fermion with wave vector k, d]^ cre- 
ates an impurity atom with wave vector k, and 6]^ creates 
a closed-channel molecule of center-of-mass wave vec- 
tor k. The Mk's obey the usual fermionic anticommu- 
tation relations, {■Uki'^k/} ~ '^kk'- Since there is at most 
one impurity atom and one closed-channel molecule in 
our system, their exchange symmetry is irrelevant. We 
can imagine for example that the impurity is a fermion, 
and that the closed-channel molecule is a boson, so 
that {dki^k/} — '^kk' and [6k, ^k'] = f^kk'- The single- 
particle kinetic energies are £k = -^^^ for the fermions, 
Ey: = for the impurity, and 



M 
m ' 



(2) 



is the impurity-to-fermion mass ratio. E^oi is the in- 
ternal energy of the closed-channel molecule. The last 
contribution in Eq. ([T]) corresponds to the coupling be- 
tween open and closed channels with a coupling constant 
A. It contains the cut-off function x depending on the 
relative wave vector for a fermion of wave vector k and 
the impurity of wave vector k': 



M 



(3) 



where fj, = mM / (m -I- M) is the reduced mass, x is a 
real and isotropic function of kj-ci. Its specific momentum 
dependence is here irrelevant, provided that x(ki.ei) tends 
to unity for fcrei 0, and tends rapidly to zero for fcroi — >■ 
+00, so as to make the Hamiltonian problem well defined. 

In practice, in all the explicit calculations to come, we 
shall take the infinite cut-off limit, that is the limit where 
the momentum width of the function x tends to infinity 
(so that x(krci) tends to unity at all fixed kioi). In this 
limit, A is fixed but -Bmoi is continuously adjusted so that 
the interaction is finally characterized by the usual two 
physical parameters, the s-wave scattering length a (of 
arbitrary sign) and the always positive Feshbach length 
R* In terms of the effective coupling constant. 
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the internal energy of the molecule is adjusted as (see e.g. 
fH): 



E, 



mol 



A2 



1 

9 



(2^ 



X'(k) 



2^1 



(5) 



and the interchannel coupling amplitude is related to the 
Feshbach length by 



AV' 



(6) 



Interestingly, i?* can be related to experimental parame- 
ters characterizing the Feshbach resonance, in particular 
its magnetic width Ai? |27j . In the present infinite cut-off 
limit, note that the s-wave scattering amplitude between 
a fermion and the impurity in the center-of-mass frame 
is simply 



/fercl ^ 1 



1 



+ ikrc\ + 



(7) 



so that the effective range is Tg = — 2i?*. By inspection 
of the poles of the scattering amplitude one readily sees 
that the model supports a two-body bound state in free 
space iff a > [28] . 

To be complete, let us recall the applicability of this 
narrow Feshbach resonance model. In an experiment, the 
direct interaction between atoms in the open channel is 
not zero, it has a finite scattering length, the so-called 
background scattering length abg (of arbitrary sign) , and 
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a finite range, the so-called van der Waals length b (of 
positive sign) . If Obg and b are much smaller in absolute 
value than the scattering length a, the Feshbach length 
i?* , and the inverse of the typical relative momentum of 
two particles in the gas, one can formally take the limit 
flbg , 6 — > for fixed values of a and i?* , and one then 
expects to reach the present narrow Feshbach resonance 
model [29|. 

III. POLARONIC AND DIMERONIC 
VARIATIONAL ANSATZ 

A. Physical motivation for the ansatz 

The goal is to determine the ground state energy of a 
single impurity interacting with A'^ fermions. In the two- 
channel model, this means that the system will in general 
coherently populate two sectors, the one with N fermions 
and the impurity, and the one with — 1 fermions and 
one closed-channel molecule. 

The exactly solvable artificial case of A = (perfectly 
decoupled open and closed channels) for a finite momen- 
tum cut-off can bring some inspiration. The ground state 
in the iV-fermion sector is then simply the Fermi sea (ab- 
breviated as "FS" ) of the N fermions with the impurity 
in the plane wave of zero wave vector. The resulting en- 
ergy is then approximately, in the thermodynamic limit: 

EpsiN) ~ InEf, (8) 
5 

where we have introduced the usual Fermi energy, ex- 
pressed in terms of the Fermi momentum fcp and of the 
atomic density p ~ N/V as 

Ef = ^ with fcp = (67r2p)i/3. (9) 

We shall use E-ps{N) as the reference energy in what 
follows. The ground state in the sector with iV — 1 
fermions and one closed-channel molecule is simply the 
Fermi sea of the — 1 fermions plus the molecule with a 
zero momentum center of mass, with a resulting energy 
Efs{N — 1) -I- Bnioi- One thus has two branches of so- 
lutions, with an energy difference in the thermodynamic 
limit 

[EpsiN - 1) + Smoi] - Efs{N) ~ E^oi - Ep. (10) 

One or the other is the ground branch, depending on the 
sign of -Emoi — Ep, which is of course artificial. 

We now turn back to the interacting case, where A > 
and the infinite momentum cut-off limit is taken. As in 
the previous A = case, there will be two branches, the 
main difference being now that the Fermi sea is depleted 
by an arbitrary number of particle-hole excitations, mak- 
ing the problem non-trivial. The first branch, emerging 
from the A^-fermion plus impurity state, was intensively 
studied in the case of zero-range single-channel models 



and is called the polaronic branch [sl-flTI. [sol Isij . The sec- 
ond branch, emerging from the A'^ — 1-fermion plus closed- 
channel- molecule state, also exists in single-channel mod- 
els flol.IT2l.l30j. It was studied in a two-channel model in 
(32i , .33| but only in the limit of a very broad and zero- 
range Feshbach resonance, where it becomes equivalent 
to the zero-range single-channel case already considered 
in [10, 12,^]. In what follows, we shall call this branch 
the dimeronic branch. 

Following these references, we explore the two branches 
on a narrow Feshbach resonance using simple variational 
ansatz with at most one pair of particle-hole excitations 
of the Fermi sea. To be able to determine the effective 
masses of the polaron and of the dimeron, we consider 
the general case of a non-zero total momentum for the 
system: 

P = hK. (11) 
For the polaron we thus take the ansatz [33 | 

l^poi(P)) - + -^q^K+q^q 
q 

+ E '^kqdl+q-k'^k^q) |FS : A^) . (12) 

k,q 

The prime above the summation symbol means that the 
sum is restricted to q belonging to the Fermi sea of A^ 
fermions, and to k not belonging to that Fermi sea. For 
K = 0, the successive terms in (|T^ correspond in that 
order to the ones generated by repeated action of the 
Hamiltonian H on the A = polaronic ground state. 
Similarly, for the dimeron we take the ansatz [sj] 

/ I 

IV'dim(P)) = ('7 4 + E^krfl_k4 + E'?kq^K+q-k4"q 
k k,q 

/ 

+ E ^k'kq4+q-k-k'4'4^q) |FS : A^ - 1), (13) 
k',k,q 

where the symbols ' above the sums indicate that they are 
restricted to k, k' out of the Fermi sea of A^ — 1 fermions, 
and to q inside that Fermi sea. Here also, for K = 0, 
the order of the terms corresponds to successive actions 
of H on the A = dimeronic ground state. Whereas 
the ansatz (|13p already appeared in [33, in the case 
P = 0, the ansatz is to our knowledge new. 

B. The integral equations to be solved 

One then has to minimize the expectation value of H 
within each ansatz, with respect to the variational pa- 
rameters that are the complex number (f) and the func- 
tions 0q and 0kq for (|12p . and that are the complex num- 
ber 7y and the functions 77k, f?kq, ??k'kq for ^61i] . As the 
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calculations for the polaron are similar but not identical 
to the ones in Q we give details of the derivation in the 
appendi x lAl For the dimeron the calculations are quite 
close to f32l. [33! so we directly give the resulting integral 
equation in the thermodynamic limit. 

For the polaronic case, one can express 0kq in terms 
of (j)q and then i/iq in terms of 0, so that one is left with 
a scalar implicit equation for the polaron energy counted 
with respect to the A^-fermion Fermi sea energy, 

Ai;poi(P) = -Bpoi(P) - EpsiN). (14) 
The implicit equation in the thermodynamic limit is 

^^'■''■''^^■'+/'(0 i>,[A^,'.,(p).pr 

where the prime on the integral over q means that it is 
restricted to the Fermi sea q < kp. The function of the 
energy appearing in the denominator of the integrand is 



and 

akciiE.P) ^ Dq{E + EF -ek,P -hk) (21) 

where the function Dq defined in Eq. (fTB)) is calculated ex- 
plicitly in the Appendix [BJ Solving the integral equation 
psp for a non-identically-zero function jykq is equivalent 
to finding the energy E ~ Ai?dim(P) for a given P such 
that an eigenvalue \{E, P) of the operator A'I[E, P] van- 
ishes: 

A[ASdi,n(P),P] =0. (22) 

Since the operator M[i?,P] is real symmetric for real 
E and P, its eigenvalues are real, and may be found 
numerically quite efficiently. The implicit equation (j22p 
is then solved numerically by dichotomy or by Newton 
method [1^1 after a Gauss-Legendre discretization as in 



Dq iE,P) = -- 



1 /ifcp fi^R* 



9 



1 



2/i 



(27r)3 V £^K+q-k' + ek' - Eq - £^ 



(16) 



where the prime on the integral over k' means that it is 
restricted to k' > kp. As shown in the Appendix [B] the 
integral in (|16l) can be calculated explicitly. 

For the dimeronic case, one can express 77, r/q and rykk'q 
in terms of Tjkq- We introduce the energy of the dimeron 
counted with respect to the energy of the Fermi sea of N 
fermions: 



AE^i^iP) = E^,„,{P) ~ EpsiN). 



(17) 



In the thermodynamic limit, this energy difference is such 
that the following integral equation for Ty^q admits a non- 
identically-zero solution: 



[AE^imiP), P; k, q; k', q']7yk'q' = 0, (18) 



(27r)e 



to be satisfied for all k such that k > k-p and for all q 
such that q < k-p. As previously mentioned, the prime 
on the integral sign in Eq.(fT8|) means that the integration 
over q' is restricted to q' < kp and the integration over 
k' is restricted to k' > kp. The energy- dependent kernel 
is given by 



X[i?,P;k,q;k',q'] = 



(27r)3j(k-k') (2^)35(q-q') 



Fkoo(£^,P) 



-(27r)6^(k-k')<5(q~q')akq(£;, P) + 



i^k'kq(£^,P) 
1 



(aooFkooi^k'oo)(£^,P) 



C. A discrete state coupled to a continuum 

Before turning to the solution of the integral equation 
it is useful to have a more global view of the polaronic 
and dimeronic problem. In the ansatz (IT^ and it 
is apparent that the first contribution is a single discrete 
state, while the other contributions form in the thermo- 
dynamic limit a continuum indexed by vectors q, k span- 
ning a continuous space. The problem then naturally 
falls into the well studied class of a discrete state coupled 
to a continuum [36| . There are then two cases. In the 
first case the discrete state is expelled from the contin- 
uum and remains a discrete state even in the presence of 
the coupling. In the second case the coupling is too weak 
to expel the discrete state which is then washed out in the 
continuum and gives rise to a resonance with a complex 
energy. This resonance can still be of physical interest 
if the imaginary part of its energy is much smaller than 
the corresponding real part. We shall encounter these 
two situations, sketched in Fig. [1] with the polaron and 
the dimeron. 

In practice one can even calculate the lower border of 
the continuum A_Ep°f*(P) for the polaron. One has to 
minimize the energy within the variational ansatz (jl2p 
setting = 0. A trivial upper bound is found by set- 
ting also 0q = and, ioi K < 2kp, we get simply 
Ai?p°f'(P) < 0. However, including the coupling be- 
tween (j}q and (/ikq, we found that the actual lower bor- 
der of the continuum can be negative for some values of 
the physical parameters kpa, kpR^, and the mass ratio 
r = M /m. The reasoning is as follows: Let us assume 
that for a given q the continuum eigenenergy AE^, is 
negative, then it is given by the implicit equation [37| 



(19) 



^q(Ai?p^ol.P)=0- 



(23) 



where 



i^k'kq(£^,P) = i^+£^F-i^K+q-k-k'-£k-ek'+eq, (20) 



Since Dq(Ai?,P) is an increasing function of AE we 
can immediately distinguish two cases. In the first case. 
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FIG. 1: The polaronic and dimeronic problems correspond to 
a discrete state coupled to a continuum, (a) When the dis- 
crete state is on the border of the continuum it is expelled 
by the coupling and remains a discrete state. (b,c) When the 
unperturbed energy of the discrete state is inside the contin- 
uum the discrete state is either (b) expelled or (c) washed 
out in the continuum where it gives rise to a resonance with 
a complex energy. The symbol + corresponds to a discrete 
state while the symbol * corresponds to a resonance. 



Z?q(0 , P) is non positive for all q within the Fermi sea. 



Therefore Eq. (p3)) cannot have a negative solution. We 
conclude that the border of the continuum is simply 
A£'p°'j''(P) = 0. In the second case, the lower border 
of the continuum A£'p°f'(P) is negative, it is obtained 
simply by minimizing the negative solutions Ai^p^j over 
q inside the Fermi sea. Let us now illustrate this formal 
discussion with the case P = 0. We then always found 
numerically that Dq(AE, 0) for a fixed AE is maximal 
on the Fermi surface [38| . The lower border of the contin- 
uum is then obtained by solving Eq. (j23|) with q = kp. By 
expressing analytically the condition Dkpe^{0~ ,0) > 0, 
we get the condition for this lower border to be negative: 



1 

kpa 



1 + r 



fcpi?* 1 - 



Inr > 0. (24) 



We now discuss the case for the dimeron setting rj = 
in Eq. (jl3p and minimizing the energy. Once again, a 
straightforward upper bound for the lower border of the 
continuum AE^°^(P) is found by setting also ?7k and rji^q 
to zero and, ioi K <3k^,we find AE^?^^{P) < 0. A less 
trivial upper bound [39| can be obtained by taking into 
account the coupling of ?7k'kq to r}k and ?7kq, where we 
realized that the particle created by uj^. is simply a spec- 
tator in the energy minimization and can be factorized 
out. This leaves us with an ansatz to minimize which is 
formally equivalent to the polaronic ansatz (I12[) with a 
fixed total wavevector K — k except for the replacement 
of |FS : N) with ul\FS : N-1). To obtain this upper 
bound it remains to minimize over k out of the Fermi 



sea 



dim (P) < , mill 



AE^l,{P-h].) + ^-E^ 



, (25) 



where AE^^^i{P) is the absolute ground state of the po- 
laronic ansatz (fT2|) (that is taken over the discrete and 
continuous spectra). We now discuss this result for van- 
ishing total momentum P = 0. In this situation we ob- 
served numerically that the global minimum of Eq. (|25|) 
is obtained at the Fermi surface so that 



AE^ZHO) < AE^l.ihkFe, 



(26) 



and can be negative depending on the physical parame- 
ters kpa, kpR^ and the mass ratio r = M/ra. Fig. [2] sum- 
marizes the discussion above, where for a given choice of 
the physical parameters we plot the lower border of the 
continuum for the polaron as well as the upper bound for 
the lower border of the P = dimeron continuum. 

To conclude we note that the qualitative features that 
we have discussed apply only for the specific variational 
ansatz ([T2|) and that we have considered. In an 
exact theory with an arbitrary number of particle-hole 
excitations the polaronic and dimeronic solutions at P 7^ 
are actually resonances rather than stationary states 
because they can radiate particle-hole pairs of arbitrary 
small energies [40|. This implies that the gap between 
the discrete state and the continuum sketched in Fig. [T] 
will close. 



IV. PROPERTIES OF THE TWO BRANCHES 
AT P = 

We calculate the energies of the two branches by solv- 
ing Eq. ([T5|) for the polaron and Eq. (1^^ for the dimeron, 
at P = 0. In both cases we assume spherical symmetry 
in the system. This amounts to assuming that the func- 
tion T^kq appearing in the integral equation (|18p has total 
angular momentum £ = 0, so that we can write 



(27) 



with being the angle between the vectors k and q. 
For the polaronic branch, assuming spherical symmetry 
means that the function 0q appearing in the ansatz (jl2p 
is spherically symmetric, as explained in the Appendix 

In Fig. 13] we plot the phase diagram of the system in 
the {k-p R*)-{'i / k-po) plane for three different values of the 
mass ratio: From top to bottom respectively r — 0.1505 
corresponding to a ^Li impurity into a Fermi sea of ^"^K, 
r = 1, and r — 6.6439, which corresponds to a ^'^K im- 
purity into a Fermi sea of ^Li. In the polaronic sector 
(left panels) we identify three different regions: (i) the 
region below the black solid line is where the polaron 
is the ground branch, (ii) the region between the black 
solid line and the lower boundary of the gray-shaded area 



6 



-9.2 r 



Polaron to dimeron transition 



lower border of the 
_9 3 _ ^\ polaron continuum 



-9.4- 



-9.5 - polaron 



-9.6- 



upper bound for the 
lower border of the K=0 
dimeron continuuin 



dimeron for K=0 



-9.7, 



0.2 0.4 0.6 O.i 



1.2 



1.4 



FIG. 2: (Color online) Polaronic energy as a function of the 
total momentum P = KK and dimeronic energy at P = 0. 
The thick black line is the solution of the implicit Eq. (|15|) 
over the interval of K where it admits a real solution, which 
corresponds to the discrete state discussed in the text. The 
dashed black line is the lower border of the continuum for 
the polaron obtained by solving Eq. (|23|l with q on the Fermi 
surface [ssl ]. The red horizontal solid line is the solution of 
the implicit Eq. (|22p for the dimeronic energy at P = 0, 
and the red horizontal dashed line corresponds to the upper 
bound (|26p for the lower border of the P = dimeron con- 
tinuum. The physical parameters are: kpR, = 1, r = 0.1505 
and l/(fcFa) — 4. Note that for this set of parameters the 
polaron has a negative effective mass at zero momentum, and 
the upper bound for the lower border of the P = dimeron 
continuum is negative. 



is where the polaron is metastable having a positive ef- 
fective mass JTip^j > 0, (iii) the region delimited by the 
gray-shaded area is where the polaron is thermodynam- 
ically unstable due the the fact that its effective mass 
is negative Tip^^j < 0. Instead, in the dimeronic branch 
(right panels) we distinguish between four different re- 
gions: (i) the region above the black solid line is where 
the dimeron is the ground branch, (ii) the region between 
the black solid line and the red dashed line (excluding the 
gray-shaded area) is where the dimeron is metastable 
featuring a positive effective mass mjjjjj > 0, (iii) the 
gray-shaded area is where the dimeronic effective mass is 
negative m'^-^^ < and indicates a thermodynamic insta- 
bility, (iv) the region below the red dashed line is where 
the dimeronic branch becomes complex, with an energy 
Aii'dim having a positive real part and a negative imag- 
inary part. In what follows we derive and discuss the 
different features of the phase diagram shown in Fig. |3l 
The main effect of having a narrow Feshbach resonance in 
the regime fcpi?* 3> 1 is that the transition from having 
the polaron as the ground state to having the dimeron 
as the ground state takes place for negative values of the 
scattering length, whereas it takes place for positive scat- 
tering lengths in the zero-range limit kpR^, — > 0. 



As discussed in Ref. tlOi] using diagrammatic Monte 
Carlo simulations, for a unity impurity-to-fermion mass 
ratio (r = 1) and for a wide Feshbach resonance (/cpi?* = 
0) the polaron to dimeron transition is characterized by 
a crossing between the two branches, the crossing point 
being located at l/{k-pa)c — 0.90. This is shown in 
Fig.|l](bl) where we plot the energies of the two branches 
calculated with Eq. (fTS)) and ([22)) at P = 0, together with 
the Monte Carlo results of Ref. [10] . On the right hand 
side of the crossing point l/(fcFa) > l/(fcFa)c the ground 
branch is dimeronic, while on the left hand side of the 
crossing point l/{kpa) < l/{k-pa)c the ground state is 
polaronic. In Fig.|4](b2) and (b3) we show how the ener- 
gies of the polaron and of the dimeron are modified by a 
narrow resonance. In Fig. 2] (a) we plot the polaronic and 
dimeronic energies for the mass ratio r — 0.1505, which 
corresponds to a ^Li impurity into a Fermi sea of '^'^K, 
while in Fig. |4] (c) we plot the polaronic and dimeronic 
energies for the mass ratio r = 6.6439, corresponding to 
a '^'^K impurity into a Fermi sea of ^Li. To summarize, in 
Fig. [5] (a) we plot the polaron-to-dimeron crossing point 
for the above mass ratios and for < kpR^, < 10. For a 
broad resonance kY■R^, =0 (a) we find that the crossing 
point is positive for the whole range of plotted mass ra- 
tios r. At r = 0.1505, for the whole range of Feshbach 
lengths we find that l/{kpa)c is positive on the range 
covered by the figure, but for larger fcpi?* (not shown) 
we find that the crossing point eventually becomes neg- 
ative. For r — 1 and r — 6.6439 the crossing point is a 
monotonically decreasing function of fcpi?* with a faster 
decrease for a higher mass ratio r. For completeness, in 
Fig. [5] (b) we plot the crossing point as a function of the 
mass ratio (r — M/m), for different values of the Fesh- 
bach length. 

Finally, let us discuss the validity of our variational 
approach where we cut the Hilbert space at one pair of 
particle-hole excitations both for the polaronic and for 
the dimeronic ansatz. For the polaronic ansatz, it was 
shown in the second paper of Ref. [lli] for a broad reso- 
nance and r — 1 that going from one pair of particle-hole 
excitations to two pairs brings typically a 10~^ correc- 
tion to the energy in units of Ep. For the dimeronic 
ansatz one may consider the zeroth order approximation, 
in which no particle-hole excitations are present, and cal- 
culate the correction induced by the first order approxi- 
mation Eq. (|13p in which one pair of particle-hole exci- 
tations is present. Setting yykq — jyk'kq = into Eq. ([T3|) 
leads to a Cooper-type ansatz given by the following ex- 
pression 



|^Cooper(P)> = (^^lc + E^k4-k4)|FS : iV - 1 



(28) 

which, after minimization of the expectation value of H, 
at P = leads to the Cooper's energy Ai?coopcr(0) = 



7 




FIG. 3: (Color online) Phase diagram of the system at P = and zero total angular momentum in the (fcF-R*)-(l/fcFa) 
plane. The left panels describe the polaronic sector while the right panels belong to the dimeronic sector. GS indicates that 
the corresponding branch is the ground state, and the black solid line marks the polaron-dimeron transition, rn'^^if^^^^-^ > 
indicates that the corresponding branch is metastable with a positive effective mass, the gray-shaded area is the region where 
the corresponding branch is thermodynamically unstable due to a negative effective mass rn^^if^^-^^-^^ < 0, the red solid line 
indicates where the inverse effective mass is zero. For the dimeron, Im(Ai5dim) < indicates that its energy is complex with 
negative imaginary part (which implies a positive real part, so that the dashed red lines correspond to A_Edim = 0). Note that 
in the ground state of our variational ansatz the effective mass is always positive, as expected, except in a narrow region in the 
upper-right panel (see subsection IIV F[) . 



-^Cooper -E-FS {N) given by the implicit equation 

aoo[AScooper(0), 0] = ( A^Coopcr (0) + Sp, 0) = 0, 

(29) 

where the definition of a in Eq. ([?T|) was used. We can 
thus estimate the accuracy of the dimeronic ansatz (jl3p 
by calculating the relative energy difference between the 
Cooper and the dimeronic ansatz at the energy crossing 
point between the polaron and the dimeron: 



S = 



|A^dim(0)- AScoopcr(0)| 



(30) 



|A£:dim(0)| /feFa=(fcFa)e 

For a Feshbach length fcpi?, = 10, and for a mass ratio 



r — M/m ranging from 1/8 to 8, we have found a max- 
imal value for S equal to 10~^, indicating that we have 
already reached a good convergence of the dimeronic en- 
ergy A£'dim(0) in terms of the number of particle-hole 
excitations. 



B. Effective mass 

The polaron and dimeron effective masses m* are de- 
fined by the expansions for a vanishing total momen- 
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M/m = 0.1505 



M/m = 1 



M/m = 6.6439 




FIG. 4: (Color online) Crossing between the polaronic energy 
Eq. p5|l and dimeronic energy Eq. (|22|) . for different values 
of the mass ratio: r — 0.1505 (column a), r = 1 (column b), 
r = 6.6439 (column c), and of the Feshbach length, respec- 
tively kpRt — (row 1), kpR, = 1 (row 2) and kpR, = 10 
(row 3). Thick lines correspond to the ground state branch 
while the thin lines indicate that the corresponding branch is 
metastable. Instead the dashed lines correspond to a nega- 
tive value of the effective mass. The red circles (dimeron) and 
black triangles (polaron) with error bars are the Monte Carlo 
results of Ref. . 



turn P: 



A£;(P) = A£;(o) 



2m* 



0(p4 



(31) 



The calculation technique for these effective mass are pre- 
sented in detail in the Appendix [Cj A negative effective 
mass indicates in our model a thermodynamic instability: 
It cannot manifest itself for an isolated evolution since 
our Hamiltonian conserves the total momentum, but it 
will manifest itself if the system is weakly coupled to a 
thermal bath, or in a harmonically trapped system, as in 
real experiments, since P is then not anymore a constant 
of motion and will slowly drift away from to reduce the 
system energy. 

In Fig. [S] we plot both the dimeronic and the polaronic 
inverse effective mass as a function of l/{kpa) for differ- 
ent values of the Feshbach length and different values of 
the mass ratio. Notice that in panel (bl), corresponding 
to a broad resonance and r = 1, we also plot the Monte 
Carlo data from Ref. ; while the dimeronic inverse ef- 
fective mass shows a very good agreement with the Monte 
Carlo data, the polaronic inverse effective mass is slightly 
larger than the corresponding Monte Carlo value. As ex- 
plained in the second paper of Ref. [11] this difference is 
due to the truncation of the Hilbert space spanned by the 
polaronic ansatz to at most one pair of particle-hole exci- 
tations. In the same reference it was shown that includ- 
ing two pairs of particle-hole excitations in the polaronic 
ansatz leads to an excellent agreement with the Monte 
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FIG. 5: Dimeron-polaron crossing point l/{kFa)c as a func- 
tion of various parameters, (a) As a function of the Feshbach 
length Rt, for three values of the mass ratio {r = 0.1505, 
r = 1, r = 6.6439 from top to bottom), (b) As a function 
of the mass ratio r = M/m for three values of the Feshbach 
length: (bl) kpR, = 0, (b2) kpR* = 1, (b3) kpR, = 10. For 
each of the figures, l/(/cFa)c is the line in the plane (kpR*)- 
(1/fcpa) or (M/m)-(l/fcpa) that separates the region where 
the dimeron is the ground state (above that line) from the re- 
gion where the polaron is the ground state (below that line). 
Solid lines: Numerical solution from the full polaronic and 
dimeronic ansatz. Dashed lines: Asymptotic analytical pre- 
diction for kpRt — >■ +00 (|59[) . The dashed line is barely dis- 
tinguishable from the solid line for r = 1 in (a) and for (b3). 
Although it does not a priori make sense, we have plotted 
this asymptotic prediction even in the regimes fcpT?, < 1 or 
— 1 < l/(fcpa)c to show that it remains quite close to the 
numerical results even in these non-asymptotic regimes. 



Carlo data. However, the inclusion of two particle-hole 
excitations is beyond the scope of this work. 

For all mass ratios considered here, in the limit where 
a — > 0+ the effective mass of the dimeron tends to the 
sum of the mass of the impurity plus the mass of a 
fermion m^^^^ — m + M. On the polaronic branch, in 
the limit where a — >■ 0^ the polaron effective mass tends 
to the mass of the impurity rn^^^ — M, since as a — >■ 0~ 
the system tends to a configuration of an ideal Fermi gas 
of A'' -I- 1 non-interacting particles. In between these two 
limits, depending on the mass ratio and on the Feshbach 
length we find different values of the effective masses as 
shown in Fig. [Bl The vertical dotted lines indicate where 
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FIG. 6: (Color online) Inverse effective mass of the dimeron 
(red lines) and of the polaron (black lines), for different val- 
ues of the Feshbach length, respectively fcpT?* = (row 1), 
kpRt = 1 (row 2), fcpi?* = 10 (row 3), and different values 
of the mass ratio: r = 0.1505 (column a), r = 1 (column b) 
and r = 6.6439 (column c). The vertical dotted lines are in 
correspondence of the crossing value l/(fcFa)c: On the left 
of this line, the ground state is polaronic while on the right 
the system has a dimeronic ground state. Thick lines corre- 
spond to the ground state, thin lines indicate where the sys- 
tem is metastable and dashed lines where the effective mass is 
negative. The red circles (dimeron) and black triangles (po- 
laron) with error bars correspond to the Monte Carlo results 
of Ref. [l^. The inverse effective mass of the polaron tends to 
m/rrip^i — >■ 1/r in the limit where a — >■ 0^, while the inverse 
effective mass of the dimeron tends to m/m*^^^ — >■ 1/(1 -I- r) 
in the limit where a — >■ 0^. The gray-shaded area in (al) in- 
dicates an instability of the system (within our P = 0, £ = 
variational calculation): The ground branch is dimeronic and 
its effective mass is negative. 
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FIG. 7: (Color online) Closed channel population (with P = 
0) as a function of l/(feFa) for different values of the mass 
ratio: r = 0.1505 (a), r = 1 (b) and r = 6.6439 (c). These 
correspond to the three bottom panels of FigU where fcp J?» = 
10. Thick lines correspond to the ground state while thin lines 
indicate that the corresponding branch is metastable. The 
vertical dotted lines are in correspondence of the polaron- 
to-dimeron crossing point: On the left the ground state is 
polaronic and on the right it is dimeronic. 



gion where the dimeron is not any more the ground state 
(on the left hand side of the vertical dotted line). How- 
ever, we observe that for the panel (al) corresponding 
to fcp-R* = 0, there exists a point \l{k^a) in the region 
where the dimeron has a lower energy than the polaron 
where the dimeron effective mass diverges, and in the re- 
gion between the divergent point ['mlm*^^^ = 0) and the 
polaron to dimeron crossing point, the dimeron effective 
mass becomes negative. We call this region the instabil- 
ity region, indicated by the gray shaded area in the panel 
(al). The question on the nature of the true ground state 
of the system in this instability region will be addressed 
in Sec. lIVF] We note however that by increasing the Fes- 
hbach length to the value /cp-R* = 1 the instability region 
disappears as can be seen in Fig. [5] (a2). 



the polaron to dimeron crossing point takes place; for 
Xjik^a) smaller than this value the ground state of the 
system is polaronic, while for l/(fcFa) larger than this 
value the ground state is dimeronic. Thick lines indicate 
that the corresponding branch is the ground state while 
thin lines represent a metastable branch. Instead, dashed 
lines are in correspondence of negative values of the effec- 
tive mass indicating an instability of the corresponding 
branch. Notice that for the panels labeled as (3), and for 
the panel (b2), the curve of the dimeronic effective mass 
has an abrupt stop at a negative value of Xjik-po). This 
point is where the dimeronic branch becomes complex 
with a negative imaginary part, as explained in subsec- 
tion [TvH 

We find that both for the polaron and for the dimeron 
there exist divergent values of the effective mass. For 
the polaronic branch, m/m*Qj = in the region where 
the polaron is not any more the ground state (on the 
right hand side of the vertical dotted line). Similarly, the 
dimeronic effective mass diverges rn,lm\^^ = in the re- 



C. Closed channel population 

An interesting quantity to consider is the mean number 
of closed channel molecules, defined as 

^cc = E(^k^k), (32) 

k 

where the expectation value is taken in the quantum state 
of the system. This quantity was recently measured in 
the BEC-BCS crossover for the s pin balanced Fermi gas 
by photoassociation techniques |41| . In our single im- 
purity problem, N^^ is less than or equal to one, and 
can also be interpreted as the occupation probability (or 
population) of the closed channel. 

When the system is prepared in an eigenstate of en- 
ergy E, and relative energy Ai? = E — Efs{N) with 
respect to the non-interacting Fermi sea, with a fixed to- 
tal momentum P if desired, a direct application of the 
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Hellmann-Feynman theorem gives [42 



1/k^a = 

F 



dAE 
dE, 



inol 



(33) 



the derivative being taken for fixed values of the other 
parameters of the Hamiltonian. Taking as a parameter 
the inverse scattering length rather than i?moi, as in [i^ . 
and using Eqs.(l5]) and (j6]) we can rewrite Ncc in the more 
convenient form 



dAE 2^iR^ 
d{l/a) ft2 



(34) 



This shows that A'cc vanishes in the broad resonance limit 
i?* — >■ 0, where N^c/R* is proportional to the so-called 
"contac t" E3II44II first introduced in one-dimensional sys- 
tems in [4^,|4^. Note that, within our variational ansatz, 



one has also Nr. 



I'yr+Ek.qi'^kqi 

-^cc = X]q l<^qP for the polaron, which are less convenient 
to evaluate than 

A straightforward calculation is to find A'cc, both for 
the polaron and for the dimeron, in the limit a — > 0+ 
but keeping R^, fixed. In this limit both the polaronic 
and the dimeronic energies are asymptotically equivalent 
to the energy of a dimer in free space E^ = —h^q^^i/2fi, 
with 



for the dimeron and 



roi given in 28]. One can also easily check that 
1/y/R^, so that in this simple limit the closed 



channel population becomes unity: 
lim iVcc = 1, 



(35) 



independently of the mass ratio as long as the Feshbach 
length is finite. Therefore, starting from this limit by 
increasing the value of the scattering length A^'cc must 
decrease and one has to numerically calculate its value 
using Eq. (|34|) . In Fig. [7] we plot the closed channel 
molecule population both for the polaron and for the 
dimeron in correspondence of the polaron to dimeron 
crossing at kpR^, — 10 and for different values of the 
mass ratio. The three panels (a), (b) and (c) correspond 
to the three lower panels of Fig. S) Notice that for the 
three cases considered here the polaronic closed chan- 
nel population is always smaller than the corresponding 
dimeronic closed channel population. Furthermore, as we 
will derive analytically in Sec.|Vl the dimeronic (resp. po- 
laronic) closed channel population tends to unity (resp. 
zero) at the crossing for kpR* — >■ -|-oo. 



D. Right at resonance 

An interesting point in the phase diagram is to look 
what happens at resonance, i.e at l/(fcFa) — 0. In Fig. [5] 
and Fig. iniwe plot respectively the energy of the polaron 
and the energy of the dimeron as a function of kpR^, 
for different values of the mass ratio. The thick lines 
indicate where the polaronic (resp. dimeronic) branch is 
the ground state. Instead the thin lines indicate that 




FIG. 8: Energy of the polaron calculated at l/(k-pa) = 0, 
as a function of kpR* for different values of the mass ra- 
tio: r = 1, r — 0.1505, and r = 6.6439. Thick lines indi- 
cate a polaronic ground state, while thin lines indicate that 
the polaronic branch is metastable and the ground state is 
dimeronic. Dashed lines indicate that the effective mass is 
negative. Vertical dotted lines indicate the crossing point 
between the dimeronic and polaronic branches (from left to 
right: r in decreasing order). The inset is a magnification. 



the corresponding branch is metastable. At fcpi?* = 
the ground state belongs to the polaronic branch for the 
three mass ratio considered here. When one increases 
fcp-R* the ground state eventually becomes dimeronic, at 
a value of /cp^* depending on the mass ratios. For r = 1, 
the transition occurs at kpR^, ~ 1.74; for r = 6.6439 it 
occurs at kpR^, ~ 0.31 (see inset of Fig.[5|); for r = 0.1505, 
it occurs at the large value fcpi?* ~ 18.09. As can be seen 
from the inset of Fig. [9l for i?, 00, the closed channel 
population saturates to unity. Also, the dimeronic energy 
A£'dim(0) tends to — i?p, as readily obtained within the 



Cooper approximation (|29|) [4 



E. Dimeronic branch: Analytical continuation 

To calculate the energy of the dimeronic branch at zero 
momentum P = 0, one has to evaluate the integral in 
Eq. jSl]) for E = Ai?dim(0). As explained in the Ap- 
pendix |B] this integral is well defined as long as E' < 0; 
even in the worst case E' = it is found to be convergent. 
However, we always found that the relative dimeronic en- 
ergy Ai5diin(0) calculated with Eq. ([22|) reaches E = 
at some value l/(/cpa)j; with a finite derivative daE > 
in that point. This is shown for example in Fig. [10] (a) 
with the full line vanishing at l/{k-pa)x — —0.35. It is 
therefore natural to ask if it is possible to analytically 
continue this curve and allow for complex values of the 
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1/k^a = 

F 




kpR* 

FIG. 9; (Color online) Energy of the dimeron calculated at 
l/(fcFffl) = 0, as a function of kpRt for different values of the 
mass ratio: r — 1, r = 0.1505, and r = 6.6439. Thick lines 
indicate a dimeronic ground state, while thin lines indicate 
that the dimeronic branch is metastable and the ground state 
is polaronic. An abrupt stop of the energy lines indicate that 
the dimeronic branch becomes complex. Vertical dotted lines 
indicate the crossing point between the dimeronic and pola- 
ronic branches (from left to right: r in decreasing order). The 
inset shows the corresponding closed channel occupations for 
the three different mass ratios; the curves overlap and cannot 
be distinguished at the scale of the figure. 



dimeronic energy: 

AEdnr.iO)=Ea-i — , where -Bo,r>0, (36) 

with F being the decay rate. Such complex energy so- 
lutions do not correspond anymore to eigenstates of the 
Hamiltonian, that is to stationary states, but they are 
called resonances. As explained in Sec. IIII C| this cor- 
responds to the case where the discrete state [the first 
component of the ansatz (fTB]) . here with K = 0] washes 
out in the continuum [the last component of the ansatz 
(I13p . parameterized by the continuous variables k, k', q]. 
If one prepares the system in the complex energy state, 
it will decay from that state into the continuum, with 
an exponential law exp(—Gammat) for the population, 
even though the system is isolated. At times ^ l/F, 
we expect that the impurity momentum will be given by 
q — k — k' as in the last component of (|13p and will have 
a broad probability distribution. 

We therefore have to solve the integral equation ^TE\\ 
with a complex energy-dependent kernel given in equa- 
tion Allowing E to be complex in the function 
Fk']iq{E,0) appearing in the kernel is straightforward, 
since apart from the pole at i? = it is continuous over 
the whole complex plane. Instead Q!kq(i5, 0) is a discon- 
tinuous function of E and has a branch cut at most given 
by the interval [0, -t-oo) on the real axis as explained in 
the Appendix [B) Therefore for all k > kp and q < k-p 



we perform an analytic continuation q;J^™'(_E, 0) through 
the branch cut from above so that 

^Hm [akq(So + ie, 0) - <q"*(So - ie, 0)] = 0, (37) 

and in the Appendix |B] we explicitly show how to satisfy 
this condition. 

We then solve at P = the implicit equation (1^^ 
for a complex Ai^dim with Newton's method 48]. In 
Fig.[Tn]we show a particular case for r = 1 and fcpi?* — 1, 
where the dimeronic energy becomes complex as soon as 
it crosses the zero threshold, which corresponds to the 
lower border of the dimeron continuum. This situation 
corresponds to the scenario of Fig. [TJ:. Notice that the 
imaginary part showed in panel (b) is about two orders 
of magnitude smaller than its corresponding real part in 
(a), which indicates a long lifetime of the dimeron. In 
addition, in the limit where a — > 0~ we find that the real 
part saturates to the Cooper's energy value E-p/r while 
the imaginary part vanishes: 

Re{A^dim(0)} Ep/r, (38) 

Im{Ai?dim(0)} 0-. (39) 

Even if not shown wc find for the same parameters many 
other complex solutions to the implicit equation (j22p 
which have a larger real part; the larger the real part, 
the larger the absolute value of the imaginary part. Fi- 
nally we have observed similar behaviors also for differ- 
ent values of the Feshbach length, in particular also for 
kpRt, = 0, and different values of the mass ratio r. 

F. Variational ansatz with ^ = 1 at P = 

In Fig.[6l for a mass ratio r = 0.1505, it was observed 
that, in a narrow region of parameters (see the shaded 
area), our P = and t — Q variational ground state is 
dimeronic with a negative effective mass and is unsta- 
ble. It can thus not be the absolute ground state of the 
system. By further decreasing the mass ratio, the insta- 
bility region becomes larger as can be seen in Fig. [11] for 
kpR^ — 0, whereas it is absent for kpR^ = 1 as in Fig. [6] 

Within the instability region, there are three possible 
scenarios: (i) the absolute ground state of the system is at 
some non-zero momentum Pq (as suggested by the neg- 
ative effective mass in P = 0), (ii) the absolute ground 
state has P = with a non-zero angular momentum, (iii) 
the absolute ground state is of a nature not captured by 
the variational ansatz (fT5|) . Exploring the scenarios (i) 
and (iii) is beyond the scope of this paper. We have 
unsuccessfully explored the scenario (ii) by solving the 
integral equation (|18p when the function yykq has a total 
angular momentum i — I (rather than i = elsewhere in 
this paper), both for even and odd parities of that func- 
tion, as detailed in the Appendix |D] For r = 1/8 and 
r = 0.1505, with fepi?* = 0, we found that, within the in- 
stability region, the £ = I dimeronic solution has a larger 
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FIG. 10: (Color online) Energy of the dimeronic branch 
A£'dim(0) = £'dim(0) - Efs{N) at P = 0, calculated for 
equal masses r — 1 and at kpR, — 1 (in red): (a) real 
part, (b) imaginary part (with error bars indicating the nu- 
merical uncertainty due to discretization of the momentum 
variables). The real negative branch for Ai5dim(0) is either 
below (red thick line) or above (red thin line) the polaron en- 
ergy Ai5poi(0) (black line). As soon as Ai5dim(0) crosses the 
zero threshold the dimeronic energy becomes complex with a 
small imaginary part that tends to 0~ when a — >■ 0~; in (a) 
the real part of the dimeronic energy is then indicated with a 
red dashed line. 



energy than the £ = one. The nature of the ground 
state of the system in this instability region remains an 
open question. 

Finally we checked the polaron to dimeron crossing 
points with the £ — 1 ansatz for the dimeron and for the 
polaron, the latter being derived in the Appendix [X] For 
all the crossings presented in Fig.|3]we found that the £ = 
1 ansatz provide a higher energy than the corresponding 
£ = ansatz. 



V. ANALYTICAL RESULTS IN A 
NON-TRIVIAL WEAKLY ATTRACTIVE LIMIT 
AT P = 

The usual weakly attractive limit considered in the the- 
ory of quantum gases simply corresponds to a — > 0^, the 
other parameters being fixed. In our impurity problem on 
a narrow Feshbach resonance, this analytically tractable 
limit shall not give very interesting information, as the 
leading order deviation from the ideal gas then depends 
on the scattering length a only, the dependence on the 
Feshbach R^, length being thus washed out. In particular, 
this can not explain the phase diagrams of Fig. [5] 

A more interesting weakly attractive limit consists in 
taking a — > 0~ while decreasing A so that R^a is con- 
stant. In the expression of the two-body scattering am- 
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FIG. 11: (Color online) Polaron and dimeron inverse effective 
mass calculated for r = 0.125 and for two different values of 
the Feshbach length: (a) fcpi?. =0 and (b) fcpii. = 1. The 
convention for the lines is as in Fig. [6] The gray-shaded area 
indicates what we call the instability region, where the ground 
state of the system within our variational ansatz is dimeronic 
with a corresponding negative effective mass. 



plitude, one then sees that the 1/a and the k^^iR* terms 
are of the same order, for a fixed relative momentum /crd 
(of the order oikp), so that the scattering amplitude /fe^^i 
tends to zero linearly in a but still keeps a dependence on 
the relative momentum. We shall now explore this non- 
trivial weakly attractive limit and show that it provides 
a quantitative physical explanation to the numerical re- 
sults. 



A. Results for the polaron 

Let us first sketch the usual weakly attractive limit 
a — > 0^ for i?* and kp fixed, where we expect for our 
impurity problem that A_Bpoi tends to zero (we are here 
at P = 0). In the integral appearing in the definition 
(fT6)) of Dq{E) we then replace E by zero, the integral 
remains finite for all values of q in the Fermi sea, even 
for 9 — >■ fcp. Then, for a 0~, one finds Dq{E) — 1/g 
which, upon insertion in the implicit equation for AiJpoi, 
reduces to the expected mean field result 



AEr, 



P9, 



(40) 



where p is the fcrmionic gas density. However all features 
due to the two-channel model have been left out in this 
trivial limit. 

In the non-trivial weakly attractive limit, we reduce 
the scattering length to 0~ and correspondingly decrease 
the interchannel coupling A to keep ai?* fixed. More 
precisely, we define for a < 



-i-aR,y/^kF, 



(41) 
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and we take the limit a —>■ 0~ for a fixed s. Assuming 
that ASpoi tends to 0~ in that limit, so that the integral 
appearing in the definition (fT6|) of the function Dq{AE) 
remains bounded 491, we have 



1 - {sq/kpf 



(42) 



As Dq has to be negative for all q < kp, this set of 
assumptions is consistent for s < 1. After evaluation 
of the integral over q in the implicit equation (jlSp for 
Aii'poi, we obtain for s < 1: 



fi^ky 1 + r kpa 1 



llni±^-l 
2s 1 - s 



(43) 



Note that, for s <C 1, this reduces to the usual mean 
field result PO]) . This calculation also allows to obtain 
the effective mass of the polaron in the non-trivial weakly 
attractive limit: In Eq. (jCll) . one replaces gDq{AEpo\) 
by its limit (|42p to obtain [50| 



M 



2fcpa 1 + 7' 



3n 



1 



+ 0((fcFa)2). (44) 



What happens to the polaronic energy when a — !■ 0^ 
for a fixed s > 1? It is reasonable to assume that AEpo\ 
has then a finite, non-zero and negative limit that de- 
pends on s. In this case 



gDqiAEpoi) 



1 



{q/k^f 



l + r AE, 



pol 



(45) 

Let us now consider the implicit equation for AEp^i'- 
Since 1/-Dq ~ one may expect that AiSpoi tends to 
zero. . . except if the integration over q leads to a diver- 
gence that compensates for the factor g, due to the van- 
ishing of Dq at the border of the Fermi surface. At this 
order of approximation, for s > 1 the discrete state en- 
ergy AEpo\ is the solution of ZJfcpe^ {AE, 0) ~ 0. It thus 
coincides with the lower border Ai?p™' of the polaron 
continuum, for P = 0, see the discussion below Eq. (P^. 
An expansion of Ai?p°"' in powers of kpa is readily ob- 
tained. To zeroth order 



AE, 



pol 



AE' 



.cont,(0) 
pol 



1 



1 



1 



(46) 



To get the next order correction, in (|45|) one now has 
to keep the additive term — fikp / {tt"^ h"^) and the integral, 
but one can approximate AE^™"' by its zeroth order ap- 
proximation P5|) in that integral. This leads to a first 
order correction 



AK 



,cont,{l) 
pol 



rEp kpa 

l+r TTS^ 



2 + J- 



1 



1 - r^/s^ 
(l + r)2 



, 1 



(47) 



where the function is defined in the Appendix iBl 

We can however be more precise and show that the 
discrete state is expelled from the continuum, as in the 



scenario of Fig. [T]d, by an exponentially small quantity. 
Expressing the energy difference as 



AEp^i - AE;if = -E 



1 



(48) 



with S > 0, it appears that the denominator 
L'q(A£'poi, 0) does not exactly vanishes in g = fcp, but 
it reaches a value ss S linearly in k-p — q, so that the 
integral in the right-hand side of Eq. (fT5|) diverges as 
gln(l/(5). Since the left-hand side of Eq. ([T5|) is close to 
Ai?p°"' , which remains finite in the non-trivial weakly at- 
tractive limit, we find that 5 vanishes exponentially fast 
for g — >■ 0~. The final result is (we recall that s > 1) [5]| : 



5 ^ A exp 

a-¥0- 



1 



{s^ - l)/{kpa) 



(49) 



where the prefactor A can be determined analytically if 
necessary. With the same reasoning applied to Eq. (jCll) , 
we conclude that, for s > 1, the inverse effective mass 
of the polaron exponentially rapidly tends to — oo in the 
non-trivial weakly attractive limit: 



M 



loi-^"- 3(l + r)5- 



(50) 



We can also proceed with the analytical evaluation of 
the closed channel population in the polaronic branch. 
For a < 0, taking as independent parameters (s,i?*) 
rather than {a,R^,), the general expression ([M)) for the 
mean number of closed channel molecules reads 



_ 771 1 + r 3 dAE 



h^kl 



ds 



(51) 



From Eas. (|43l46l) we then find the simple results: 

^ " for s < 1, (52) 
for s > 1. (53) 



iVP°' for 

o-s-O- 



Finally, we physically explain why the results for the 
polaron depend on the location of s with respect to unity, 
using the picture of a discrete state coupled to a contin- 
uum exposed in subsection IIIICl According to Eq. (1^ , 
expressed in the non-trivial weakly attractive limit, we 
immediately find that the lower border of the polaron 
continuum is zero for s < 1 and is negative for s > 1. So 
the physical situation for s < 1 corresponds to Fig. [T^, 
and for s > 1 it corresponds to Fig. [T}d. 



B. Results for the dimeron in the Cooper 
approximation 

The same non-trivial weakly attractive limit as in the 
previous subsection may be applied to the dimeronic case. 
As shown by the numerics, the Cooper approximation is 
accurate in the relevant regime, in practice close to the 
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crossing point between the polaronic and the dimeronic 
branches. We shall provide an analytical justification 
of this fact in subsection IV CI For P = 0, the implicit 
equation for the dimeron (that we now call cooperon since 
this is in the Cooper approximation) is 



i?o(Ai;Coopcr+£^F,0) =0. 



(54) 



As long as Ai?coopcr < Ep /r, the integral appearing in 
the expression of Dq is well defined. Multiplying Eq. ([5l)) 
by g and taking the limit g — )• 0^ for a fixed s, one gets 



Cooper 



Ef 



= 0, 



which leads to the prediction 



AScooper -t A4°ip„. ^ Ef 



r 1 
1 + r 



(55) 



(56) 



This is meaningful for AE^^^^^^^. < Ep/r that is for 
s > r/(r + 1) [otherwise Ai^coopcr — > Ep/r exponen- 
tially fast]. It is then found that, in this weakly attrac- 
tive limit, the polaronic and the dimeronic energy curves 



cross when s = [r/(l -|- r)]^/^, at a vanishing energy in 
that limit, which leads to the prediction 



1 

k-pa 



1 



-kvR^ 



(57) 



This reproduces quite well the slopes of the numerical 
results of Fig. [5] To get a prediction for the intercept, one 
simply has to calculate AE'cooper to next order in kpa. In 
(IMl) one now has to keep the additive term — fik-p / {ir'^f?) 
and the integral, but one can approximate Ai^cooper by 
its zeroth order approximation (j56p in that integral. This 
leads to a first order correction 



AK 



(1) 

Cooper 



h^kl 



kpa 1 



1 + r TT 

r 



, l+r 



1 - 



In- 



By equating AS^^j^p^^. + AE^^t^^^^^. to the leading order 
expression Eq. (j43p for the polaron energy we obtain for 
the crossing point in the weakly attractive limit [52]: 



2(H-r)s 



(1) 



(58) 



1 

kpa 



r 2 , 

= fcFi?* + - < 1 

J?..^+oo l + r TT 



l + r 



5/2 



l + r 



1/2- 



In- 



1/2 



l+r 



1/2 



-0 



(59) 



This asymptotic prediction is in very good agreement 
with the numerical results, already for moderately large 
values of fcp-R*, see Fig. [51 For a mass ratio M/m = 1, the 
asymptotic prediction for fcpi?* =0 exceeds the full nu- 
merical calculation by ~ 3% only, which is of course very 
fortuitous, if one realizes that the corresponding value 
of {kpa)c is positive while the asymptotic result was ob- 
tained in a limit where a -> 0^. . .Note that the inter- 
cept, that is the constant term in Eq. ([59)) . is positive for 
all mass ratios, it diverges as 2/(37rr) for r — ^ and it 
vanishes as 3[ln(4r) — 4/3]/(7rr) for r +oo. 

The closed channel population can also be accessed for 
the dimeronic branch in the Cooper approximation, from 
the general formula ((5T|) . For s > r/(l + r) one finds that 
it tends to unity in the weakly attractive limit: 



^Cooper ^gQ^ 



C. Weakly attractive limit for the full dimeronic 
ansatz 

We now argue, in the previously considered weakly 
attractive limit (a — > 0^ with ai?* fixed), and for 
P = 0, that the energy Ai?dim corresponding to the full 
dimeronic ansatz only weakly differs from the Cooper re- 
sult Ai?coopcr, in the sense that, for fixed ai?*, 

5E = AEdin,-AEcoopcr = 0{a^). (61) 

As a consequence, the expansions (I56I58P still apply for 
Ai?dimi a-iid the asymptotic result for the crossing point 
([59| is unaffected at the considered order. For simplicity, 
we restrict to the regime where Ai^coopcr < 0, so that 
Ai?dim really exists as an energy minimum (and not as 
a complex pole in an analytic continuation, see subsec- 
tion |TVE|. This regime contains in particular the cross- 
ing point of the dimeronic and polaronic branches, as we 
have seen in subsection IV Bl 

The idea is to perform a perturbative expansion, cal- 
culating 5E to leading order in a, for a — > 0~ with s of 
Ea. (HT|) fixed. Multiplying the integral equation ([T5)) by 
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gaoo{E) / a]^q{E) for E equal to the dimeronic relative energy, we get the linear system: 
I 



5Q!Oo(Ai?dim) ?7kq 



i] 



?7k'q' 



(27r)3 (27r)3 F^oo{^Edi^)F^,oo{^Edim) 



g'^aoo{AEdim) 

m J 



'd^q' 



(2^)3 Fkoo(A^dim) J (27r)3 Fk'kq(A£;, 



'd^'k' 



??k'q 



(62) 



r 



At the energy Ai?coopcrj go^oo vanishes by definition, so 
we linearize in SE, calculating the derivative of aoo with 
respect to the energy to leading order in the weakly at- 
tractive limit: 

gaoo ( A£;coopcr + SE) ~ - ^ 6E. (63) 

The left-hand side of (p^ is thus of order ((5i?/-EF)'7kq- 
In the right-hand side of (j62p . it appears that t^akq has 
a finite limit for g 0^ with fixed s [ssj : 



ffakq(A£;Coopor) 



1 + r 
r 'E^ 



£k — £q 



^q-k 



1 + r 



(64) 



Note that this limit is non-negative, and vanishes only for 
q and k having coinciding values on the Fermi surface. 
Also the functions F in the denominators of (l62t do not 
depend on g. Taking into account the linearization 
it appears that the first term in the right-hand side of 
([S^ . which is not proportional to SE, is of order g'^. On 
the contrary, the second term in that right-hand side is 
of order gSE, with g — >■ 0~, and is negligible as compared 
to the left-hand side of (p^ . and thus as compared to the 
first term of the right-hand side. 

To leading order in g, we finally have the eigenvalue 
problem 



SE 



??kq = 



Ek — £c 



£q-k 



' d3g' d^k' 



1 

?/k'q' 



(27r)3 (2^)3 Fkoo(A4"ip„)i^k'oo(A£;, 



(0) 

Cooper'' 



(65) 



where the zeroth-order approximation for the Cooper en- 
ergy is given by ()56p . This eigenvalue problem is solved 
by integrating over k and q after division by i^koo, which 
leads to the explicit expression for the leading deviation 
of the dimeron energy from the cooperon energy: 



As promised, this vanishes as g^ in the non-trivial weakly 
attractive limit. Obviously, one also gets the correspond- 
ing leading order approximation for r?kq- Note that our 
expression for SE is negative, which was expected from 
the usual argument that the variational space for the 
dimeronic ansatz contains the cooperonic ansatz [s^ ]. In 
Fig- EH we successfully compare the energy difference 
with the corresponding energy difference calculated nu- 
merically with the variational ansatz (|13p and (f28| . 

To complete this analytical discussion of the dimeronic 
ansatz, we use the discussion of subsection IIII CI to have 
a prediction on the lower border of the dimeron contin- 
uum at P = in the non-trivial weakly attractive limit. 
Evaluating the lower border AEp°f^{h'K) of the polaron 
continuum for an arbitrary total momentum fiK can be 
done by solving Eq. (|23p in that limit and minimizing 
the root over q inside the Fermi sea: One obtains either 
KK) = for small s, or 



AiS^of 



Ai?-f(K) 



1 



-1) 



1 



1 



[{K/kpf - 2K/kY 



(67) 



when s is large enough for the right-hand side of (|67p to 
be negative. From Eq. (f26| we conclude that 



A£;™i;*(0)<min 



1 



1 + r s' 



- 1 Sf. 



(68) 



The lower border of the dimeron continuum is thus 
negative in the non-trivial weakly attractive limit for 
s > [r/(l -|-r)]^/^, which is the same range of s leading to 
a negative cooperon energy, see Eq. (|56p . For s slightly 
above [r/{\ + r)]^/^, the lowering of the dimeron con- 
tinuum border pushes down the dimeron discrete state 
energy, according to the scenario of Fig. [TJ), leading to 
its crossing with the polaron energy (that remains van- 
ishingly small for s < 1 when a — > 0~). 



VI. CONCLUSION 



SE~- 



gE^ 



' d^q d^k 
(27r)3 (27r)3 

r Ef 
1 + r s-^ 



£k — £q 



2fi 



(66) 



We have studied in detail the problem of a single impu- 
rity of mass M interacting with a spatially homogeneous 
Fermi sea of same-spin-state fermions of mass to on a s- 
wave narrow Feshbach resonance. As we have discussed, 
this problem belongs to the general class of a discrete 
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FIG. 12: (Color online) Dimeron-to-Cooperon energy differ- 
ence (|6ip divided by (fcpa)^ as a function of kpa, calculated 
for a mass ratio r = 6.6439 at s = Sc = ['"/(I + f)]^^^ where 
s is defined by Eq. (|41l) . To leading order in the non-trivial 
weakly attractive limit, s — Sc corresponds to the polaron- 
to-dimeron crossing point. The horizontal line is the analytic 
Eq. (|66p . which is reached by the full numerical calculation 
only in the limit a — >■ 0~ (with ai?* fixed). The dashed line 
is the slope of the linear terms in the numerical data. 



state coupled to a continuum. Due to the narrowness 
of the Feshbach resonance there is, in addition to the s- 
wave scattering length, a second parameter characteriz- 
ing the interaction, namely the Feshbach length i?* . Us- 
ing simple variational ansatz limited to at most one pair 
of particle-hole excitations of the Fermi sea with Fermi 
momentum kp , we have determined for the polaronic and 
dimeronic branches the phase diagram between absolute 
ground state, local minimum, thermodynamically unsta- 
ble regions (with negative effective mass), and regions 
of complex energies (with negative imaginary part). In 
particular, we have numerically calculated the polaron- 
to-dimeron crossing point l/{kpa)c, where a is the s-wave 
scattering length, as a function of i?* and for different val- 
ues of the mass ratio r — M/m. One of our main results 
is that the polaron-to-dimeron crossing point is shifted 
in the negative a region for large fcpi?*. A physical ex- 
planation of this fact is given below. Putting forward 
a non-trivial weakly attractive limit where a — >■ 0~ and 
oo, with the product ai?* being constant, we could 
obtain analytical results for the polaronic and dimeronic 
branches and finally get the asymptotic expansion (j59p 
for the crossing point. We have calculated experimentally 
accessible properties of the system, such as the effective 
mass both for the polaron and for the dimeron m*^^^, 
as well as the closed-channel molecule population A^cc for 
the two branches. 

For M/m = 0.1505 we have found a small region of pa- 
rameters where the ground state of the system restricted 
to zero total linear and angular momenta is dimeronic 



with a negative effective mass m^j^ < 0, therefore in- 
dicating an instability of the system. Extension of the 
dimeronic variational calculation to the angular momen- 
tum £ = 1 (with still P = 0) has not allowed to clarify 
the nature of the absolute ground state of the system in 
this instability region. 

We now propose a physical interpretation of the shift 
towards negative values of a of the polaron-to-dinieron 
crossing point for a narrow Feshbach resonance. First, 
in a broad Feshbach resonance, described with a single 
channel model, a two-body bound state in free space can 
exist only if a > 0. A common wisdom is that the pres- 
ence of the Fermi sea makes it more difficult for the im- 
purity to form a two-body bound state with one of the 
fermions. Such a bound state at zero momentum is in- 
deed a spatially localized object formed by the impurity 
and a fermion by pairing of their momenta ±ftk over some 
interval of values of k. The presence of the Fermi sea lim- 
its the momentum available for the considered fermion to 
be larger than hkp. 

In a narrow Feshbach resonance, described with a two- 
channel model, there is from the start a molecular state 
in the closed channel and the issue is to understand in- 
tuitively why it does not immediately give rise to a two- 
body bound state between the impurity and a fermion. 
In short, due to the interchannel coupling A, the molecu- 
lar state can dissociate into a continuum of open channel 
states where the impurity and the fermion have oppo- 
site momenta, if the molecular energy is above the dis- 
sociation threshold of the open channel. A naive appli- 
cation of this argument in free space, would lead to the 
wrong conclusion that there is a two-body bound state iff 
Emoi < 0. But the closed channel molecular state experi- 
ences a strongly negative "Lamb shift" due the coupling 
A with the vacuum fluctuation in the open channel [ssj . 
so that i?nioi is renormalized into: 



E, 



mol 



Err 



(fik ^ 2n 
X (k) 



(27r) 



h?k-^' 



(69) 



According to Eq. ([5|), -Bmoi has the simple expression 
^-moi = ^A'^/g. The condition that there is a two-body 
bound state in free space iff Erao\ < gives the correct 
condition that a > 0. In the presence of the Fermi sea 
in the open channel, there are two changes: (i) the effec- 
tive dissociation threshold is now the Fermi energy E-p 
[11], and (ii) in the "Lamb shift" ([M)) one has to re- 
strict the integration over k to /c > fcp, which makes that 
shift slightly less negative. Therefore there is a two-body 
bound state iff the resulting i^moi is less than E-p , that is 



1 2 

kpa TT m 



M 



M 



-k-pR^ 



(70) 



In the right-hand side of that equation, the first term 
originates from the modification of the "Lamb shift" by 
the Fermi sea, and the second one from the rise of the 
dissociation threshold by the Fermi energy. Remarkably, 
this simple reasoning perfectly agrees with the non-trivial 
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weakly attractive limit lowest order result (|57p . where 
kpR^ — ^ oo [13] ■ In the opposite limit of a broad Fes- 
hbach resonance, it reproduces the expected result that 
the two-body bound state cannot form unless l/k-pa is 
positive and large enough j58j . 
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Note 

While coinpleting this work we became aware of related 
studies I59ll60ll. 



Appendix A: Implicit equation for the polaron 
energy 

Here we derive the implicit equation (1151) for the po- 
laron energy A£;poi(P) = £^poi(P) - EpsiN). Equat- 
ing to zero the first order derivatives of (V'poi(P)l^^ — 
i?poi(P)|'0poi(P)) with respect to the complex conjugates 
of the variational parameters leads to the following set of 
coupled equations 



[A£;poi(P) - Ek] 



V ^ qSFSjv 



[Ai;pol(P)+£q-£k-£;K+q-k]0kq = X [k- ^ (K + q)] ,/)q 

(A2) 



A£:poi(P) + eq 

K 

M 



£K+q 

1 + r 



E„ 



A_ 

V 



L m 



)]'^+ E x[k-^(K-t-q)]</>kq}. 

(A3) 



k^FSw 



This can be reduced to a single equation for the varia- 
tional parameter (/)q: From the first equation we elimi- 
nate (j) a-nd from the second equation we eliminate 0kq, 
using the fact that the energy differences cannot vanish 
for AE'poi < 0. We then divide the third equation by 
A^ and we replace the internal energy of the molecule 
by its expression Eq.([5]). We take the thermodynamic 
limit, replacing the discrete sums over q and k by in- 
tegrals. We split the integral in Eq.(IS]) in two bits, an 
integral for k < k-p and an integral for k > kp. The bit 



with k > k-p can be collected with the integral over k 
resulting from Ea. (|A3p : This allows to directly take the 
infinite momentum cut-off limit, where x — > 1, without 
activating a divergence. Using Eq.® we finally obtain 



i?q[Ai?pol(P),P](/., 



1 fd^q' , 

AEpoi{P)~EkJ (27r)3'^'»'' 

(A4) 

where Dq[i?(P),P] is given in Eg. dTB)) . A non-trivial 



solution of Ea. (IA4[) may be easily found by making the 
change of variable 



{2n) 



3rq- 



(A5) 



Then, if i:iq[A£:poi(P), P] ^ for aU q < kp we can 



divide Eq. (|A4p by Dq and integrate on q over the Fermi 
sphere. Dividing the resulting equation by (f under the 
assumption if ^ 0,we recover the scalar implicit equation 
for the polaron energy (fT5|) . For P = the corresponding 



solution for (j)q has zero total angular momentum ^ = 
since Z3q is then rotationally invariant. Note that the 
assumption if ^ necessarily fails if one rather looks for 



(A6) 



a solution (pq with ^ > 0. Then Eq. (jA4|) reduces to 

i?q[ASpol(0),0] =0, 



for some q inside the Fermi sea. As expected, and as 
checked numerically, the solution to Ea. (jA6p has always 
a higher energy than the solution with £ = (63 | . 



Appendix B: Explicit form of the integrals in Dq, 
Qkq and in a'^^ 



First we calculate the integral 



Iq{E,P)= J 



' (fk' 
(2^ 



1 



-E^K+q-k' + Ek' — 



E 



2/i 

(Bl) 

function 

£'q(i?, P). This is equivalent to calculating the function 



that appears in the definition (|16p of the 



akq(-E, P) according to (pij) . and the reciprocal property 
holds since Dq{E, P) = aoqiE — Ep, P). To be sure that 
the integral is well defined, one has to check that the 
energy denominator in Eq. ljBip remains positive for all 
values of K (over the whole Fourier space), of q (such 
that q < kp) and of k (such that k > kp). We thus 
impose the following constraint on the energy E: 



E < inf ^K+q-k' 

K,q,k ^ 



£k' - eq = 0, 



(B2) 



where the prime on the infimum symbol recalls the con- 
straints q < kp and k > kp. Note that the condition 
i? < is not restrictive since E stands in practice for 
energy differences that are negative in the ground state 
of the system. 

The expansion of -EK+q-k' in powers of k' reveals that 
it is convenient to use spherical coordinates [k', 9' ,((>') of 
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polar axis given by the direction of K+q. The integration 
over the azimuthal angle (j)' is straightforward, and after 
the change of variables x' = cos 6' and k' = fc'/fcp we are 
left with 



to a gives 



c^+b 1 



d.Ha, b,c) = -- + -) In 



dx' 



2/ifcF 

(2^ 



+ 00 



arctan ■ 



1 



+ --t-- ^ - TT Uf 46 > 



-k'2 -v4q(P)a;'K' + Bq(£;,P) 

where we have introduced 

2^|K + q| 



^q(P) 



M 
2fi 



BqiE,P) ~ -^^[E-K+q- Sq- E]. 



(B3) 

(B4) 
(B5) 



(BIO) 



The derivative of (jB7| with respect to b gives 

a t:/ u \ 1 1 c2 + ac + 5 , 
ObJ- (a, o,c) ~ — in ■ 



arctan 



2a c? — ac + b 

^£±2_ _L arctan ?.r°., - tt ] if 46 > 



After integration over x' this leads to integrals of the 
type / dK'«;'ln(K'^ ± an' + 6), which have to be solved 
separately for the case 46 > a^ and 46 < a^ . The final 
result is given by 

I^{E, P) = ^^:F [Aq{V), Bq{E, P), 1] , (B6) 

where we have introduced the function that vanishes for 

c — > +oo: 



1 2\/ 



+1 ( In ^^±^+^^^ + In 1--+^^^ ) if46<a2. 

a2-4b\ 2c+a — Va -46 2c— a-\/a — 46/ 

(Bll) 



J^(a,6, c 

f V4b-a^ Z' 
2 I 



2c^ 



4a 



■In- 



ac 



arctan 



2c+a 
V4fc-a2 



arctan 



c"" — ac + 6 

- tt) if 46 > a^ 



In the third part of this appendix, we allow the en- 
ergy to be complex E E C and analytically continue the 
integral Iq{E,P) defined in (jB6|) . Taking the standard 
branch cut for the logarithm and square root functions 
of a complex variable, i.e. the interval (—00, 0] on the 
real axis, then the branch cut of Iq{E, P) is given by the 
interval [0,-|-oo) on the real axis. We continue Iq{E,P) 
from above over this branch cut: Given a, c e R this is 
equivalent to analytically continue J-{a, 6, c) as a function 
of 6 e C so that 

lim \T(a, bo - ie, c) - Fia, 60 + ie, c)] = 0, (B12) 

£^0+ 



Va^-4b 
4 



2c+a— x/o?— 4b 



In 



2c-a+\/a^— 4b 



if 46 < a^ 

(B7) 



In the second part of this appendix, we calculate the 
function a[^^{E,P) defined in (|C17|) . Since it is simply 
the derivative of akq(-E, P) with respect to Pz, one is 
reduced, by virtue of relation (|2T]) . to the derivative of 
£'q(i?,P) with respect to P^: 



dp^Dq{E,P) 



Ai^i?, Kz + Qz 



dpJqiE,P). (B8) 



for aU 60 = Re(6) e JR. To satisfy the condition (IBT2)) 
one has to calculate the same limits on each function 
appearing in J- but in practice there are only two types 
of functions: The first one is \n{c^ ~ ac + b) which can be 
discontinuous when Kc{c^ — ac + b) < 0: 

lim ln(c2 — ac + 60 + ie) = In jc^ — ac + 6o| ± jTr. (B13) 

We then substitute \n{c^ — ac + b) with the continuous 
function 



nh^ m + M 

From the chain rule for derivatives applied to (jB6[) : 



fia,b,c) = 

ln(c2 — ac + 6) 



if Re(c2 - ac + 6) > 



dpJq{E,P) 

( kp 



|K + q 



^ ^i^{Kz + qz) 

daJ'+dhT)[Aq{P),Bq{E,P),l]. (B9) 



Note that dbJ^ directly gives access to dEDq{E,P), and 
thus to 9_Bakq(-£': P): which is useful for a Newton search 
of the polaronic and the dimeronic energies. An explicit 
calculation of the partial derivative of (|B7|) with respect 



ln(-c2 + ac - 6) - ZTT if Rc{c^ - ac + b) <0, 

(B14) 

which satisfies /(a, 69 — ie,c) — /(a, 69 + ie,c) when 
e ^ 0+. Notice that we have conventionally chosen 
Eq. (jB14p such as to concord with the minus sign limit 
of Eq. (jB13l) in the relevant region. The second one 
is In (2c — a ± \/ a'^ — 46) which in turns has two cases. 
When Ke{a^ — 46) > 0, similarly to the previous case 
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we substitute In (2c — a ± y/a^ — 4oj with the continuous 
function 

h±{a, 6, c) = 

fin (2c-a±Va2-4fe) if Re (2c-a± Va2-46) >0 



[in (a-2cTVa2-46)±i7r if Re (2c-a± Va2-46) <0, 

(B15) 

which for e — !■ 0^ satisfies h{a,bQ — ie,c) — h{a,bQ + 
ie, c). Instead when Re(a^ — 4&) < we simply substitute 
— 46 with i^/4b — and no other arrangement has 
to be done in the logarithmic functions. 



Appendix C: Calculation of the polaron and 
dimeron effective masses 



To calculate analytically the polaron effective mass, we 
simply expand the implicit equation (jl5p up to second 
order in P = ?iK, which requires in particular an expan- 
sion of the function Z)q[Ai?poi(P), P] up to second order 
in P. After integration over the direction of q in ([T5)l, 
terms that are linear in K vanish for symmetry reasons 
and, after division by K^, one is left with 




^i^R, f'd^k' 1 



TTh^ 7 (27r)3 ^ 
(27r)3 3APF^, 



k'q _ 

(CI) 

where we have introduced the short-hand notation Vq = 
Dq[AEpa\{0), 0], and we have defined the scalar quantity 
Fk'q = i?q-k' + £k' — £q — Ai?poi(0) and the vectorial 
quantity (parallel to q for symmetry reasons): 



^^i?* q 
7r(m -f- M)h'^ 



' d^k' fi^jk' -q)/M 



^k'q 



(C2) 



One possibility to calculate mj^im ^'-'^ dimeron is 
to solve Eq. ip^ for a few values of P S [0,Pniax] and 
then extrapolate m^j^^ from Eq. ([3T|) . However, this 
procedure is computationally very costly and is not safe 
from numerical errors, in particular due to the arbi- 
trary choice of Pmax- Instead, we use a perturbative 
approach in which we take the limit P in the oper- 
ator M[_E(P),P]. Thus, a second order expansion gives 



M[£:(P),P] ~ M[E{0),0] 



p2 



2m* 



-dEM[E{0),{)] 



+ P.dp,M[E{Q), G\ + \Y1 P^Pj9p,dp^ M[^(0), 0], 



(C3) 



where the indices i,j run over the three orthogonal com- 
ponents {Px, Py, Pz) of P, and we consider the operators 
in (jC3| proportional to P to be small perturbations of 
the operator M[E{0),0]. At P 0, the unperturbed 
eigenvalue equation is given by 

M[EiO),0]\vn) = A„[£;(0),0]|i;„), (C4) 

and the shift to A„[i?(0), 0] induced by the perturbations 
is then 



A„[i;(P),P] ^ A„[i?(0),0]-h^<5Ai^)[i?(0),P], (C5) 



where each contribution is calculated using first and/or 
second order perturbation theory as follows 



a«[i?(o),p] 



p2 



{vn\dEM[EiO),0]\vn), (C6) 



JAi2)[£;(0),P] - («„|^P,9p,Af[S(0),0]|t;„) 

i 



ra^n 



A„[£;(0),0] - A,„[£;(0),0] 



, (C7) 



5\^^\e{Q),V] = M\Y.P,P,dp^dp^M[E{Q)Mvn). 

hi 

(C8) 

At P = 0, the dimeronic energy E{0) is found by 
setting to zero the minimal eigenvalue of the operator 
Af[£;(0),0], namely 



Ao[^(0),0] -0, 



(C9) 



which implies Ao[-E(P),P] = for a finite P, thus the 
effective mass is readily obtained from the implicit equa- 
tion 

^a(^'[£;(o),p] = o. (CIO) 

Specifically, for the eigenvector \vq) we use an ansatz of 
zero total angular momentum 



(k,q|z;o) = /(fc,g,0). 



(Cll) 



where 6 is the angle between the vector k and q. Without 
loosing in generality we can assume the vector P to be 
oriented along the axis, i.e. P = Pe^, and taking 
into account the symmetry of the ansatz (jClip the effect 
of the perturbations on Ao[£'(0),0] is then reduced to 
calculate 



<)|B(0).P1. 



{vo\dEM[E{0),0]\vo) (C12) 
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\{vo\dpJI[E{0),0]\v„,)\' 



(C13) 



SX^,'\e{0),P] ^ ^{vo\-^M[E{0),0]\vo). (C14) 



The effective mass of the dimeron is then easily obtained 
from Eq. (jC10|) . and we have 



\{vo\dpJI[E{0),0]\v„ 



X„,[E{0),0] 



- ("0 \-^M[E{0),0] \vo)}/{vo\dEM[E{0), 0] |t-o). 

(C15) 

While the operator M^^) [£:(0), 0] = dEM[E{0),0] is 
of straightforward calculation, the terms in bracket are 
more complicated and in what follows we explain the de- 
tailed calculation. 

First term in (|C15|) — The kernel of the operator 

is given by the expression 

A^(2)[£;,k,q,k',q'] - -(27r)6(5(k-k')<5(q-q')«kq(^'0) 
h f (2^)3^(k-kO , , (27r)3^(q-q0 , _uy, 
m\ FUiE^O) Fi,^^iE,0) 



-^^kOO ^k'OO 

where 



(aooi^kooFk'oo)"'](S,0)}, (C16) 



a;,q(£;,p) = ap,akq(s,p) 



(C17) 



can be evaluated with the formulas of appendix |B] 

Using the completeness relation 1 = the 
first term in the right hand side of Eq. (jClSP may be 
written as 



E 



\{vo\M(^HE(0),0]\ym)\' 



X„,[E{0),0] 
{vo\M^^\EiO),0] 



1 



M[E{0),0 



■M^^\EiO),0]\vo) 



(C18) 



It is not difficult to check that |a;[i?(0)]) = 
M^^^[E{0),0]\vo) may be written as 

{k,q\x[E{0)]) = ^u,[k,q,9,E{0)] + ^v,[k,q,0, EiO)]. 
k q 

(C19) 

Given that 6 is the angle between k and q, Eq. (IC19I) is 
evide ntly an odd function of (k, q) , and as explained in 
Ref. [25[ has total angular momentum i = 1. Using the 



ansatz (jClip . the functions appearing in Eq. (IC19|) are 
then given by: 



Ux{k,q,9,E) 



k 



k 



M I 47r2 _ 

1-1 (aooi^^oo^k'oo)(i^,0) 



dk'k' / dq'q Idx 

F "'0 J-i 



J-i ^ kOO 
fik'^q'^x') 



47r2 



dk'k'^i dx 



kOO^ 

f{k',q,x') 
[5g(i?)-62]3/2 



vr 



and respectively 

V,{k,q,e,E) = 



bik' -bo{E)k ^ ~a'^^{E,0), (C20) 



h 1 
M4^ 



„I2 



boiE){q-k'x')- '^L_^ bik'x 



(C21) 



where as usual x — cos 9. Where appropriate we have 
substituted 



Fk,kqiE,0) = boiE) + bi cos(^, 



(C22) 



where 



2^^ 



E + Ey 
^(k + k')q--A:fcW 



h'^q^ / 1 1 
2 Vm ^ M 



and respectively 



bi = 



M 



kk'\/l — x^\/l — x'^ , 



{C23) 



(C24) 



and we have performed the integration in cj) [63l | . 

Since the operator M[E{0),0] has a vanishing eigen- 
value Ao[i?(0),0] = it is not possible to invert it, but 
we can overcome the problem by solving 



M[EiO),0]\y)^\x[EiO)]). 



(C25) 



Notice that since the source term |x[i?(0)]) has odd parity 
and total angular momentum ^ = 1, it is easy to check 
with (|19p that also \y) must have odd parity and unit 
angular momentum. Following Ref. [l^ , this reduces the 
unknown \y) to be of the form 

(k, q|y) = ^Uy{k, q, 9) + ^«,(fc, q, 9), (C26) 
fC q 

where 9 is the angle between the vectors k and q. Writing 
this ansatz into Eq. (|C25P leads to an implicit equation 
for the energy E{0), in the unknown functions Uy{k, q, 9) 
and Vy{k, q, 9): 



k. 



{Uk + Vk-u,}[k,q,9,E{0)] 

+ ^ {[/ + 1/ y^} [fc, 0^ E{0)] = 0, (C27) 
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where 



given by 



Ukik,q,e,E) 

+ ^ / dk'k- 



j^, Uy{k,q',x') 



1 

1 

1 v/62(^y^yr^ 



a; 



/2 



Uy{k',q,x') 



Vk{k,q,0,E) = — 



ai,q{E,0)uyik,q,e), (C28) 



'„'2 



-1 Fkoo{E, 0) 



(C29) 



{27T)^S{k - k') 2ft2 fc2 1 1 



(27r)3^(q-q') 



3M i^koo(i^,0) 
2ft2 (q - k 



M 



M 



a, 



00 



1 



aoo M V Fkoo 
2^2 k • k' 

X 



k'OO 



fc2 



2fi2 



-^k'OO 



M2 FkooFk^oo 
where 



(aooFkoo^^k'oo)"'}(S,0) (C33) 



U,ik,q,9,E) 



X 



dk'k'- 



1 

1 

dx'- 



dk'k'^ 

F 

2+ 



dx'x' 



rUj,(fc',g,x'), 

(C30) 



1 f f 

V„{k,q,e,E) ^ — dk'k'^ dx 



-akci{E,O)vy{k,q,0). (C31) 

Once again x = cos 6, where appropriate we have made 
the substitution_ (IC22p and performed the azimuthal in- 
tegration in (j) [63|. 

Therefore, for a given energy E = E{0) we find the 
unknown functions Uy{k,q,0) and Vy{k,q,6) by setting 
to zero the two hues of Eg. (|C27|) separately since they 
are linearly independent [63|. Then the scalar product 
(jClSp is straightforward and is given by the integral 



[Em\y) = 



247r4 



dkk j dqq j dx 



{< [fc, q, 9, E{0)]uy{k, q, 6) + vl [k, q, 6, E^Vyik, q, 0) 



a'i^iE,P)=J2^aaE,P) = 



1 r' d^k' r 1 



2^2 (K + q - k - k')2 



t,^^{E,P) 3M Fl,^^{E,P) 



(C34) 



As previously mentioned, the prime on the integral sign 
means that the integral is restricted to q < kp and/or to 
k > kp. Eq. (|C34P may be calculated from Eq. ([21]) by 
using the simple relation 



a'LiE,P) = 



7rft4 Km 



1 



M m + My 

1 2 
dEakq{E,P) - -^^di/MdE[akqiE,P)]R,=o. 



M 



(C35) 



Appendix D: Dimeronic ansatz with £ = 1 

Here we solve the integral equation (|18p for 7/kq having 
total angular momentum i = 1, both for even and odd 
parity. Following Ref. [1^, for an even function of total 
angular momentum £ = 1 we write the following ansatz 



?7kq 



k A q 



f{k,q,0), 



(Dl) 



where is a unit vector along the z direction, A indicates 

ul[k,q,0,E{O)]vy{k,q,e)+v*[k,q,0,E{O)]uy{k,q,e)] }. vectorial product, 6 is the angle between the vectors 

-! J k and q. We insert the ansatz (|Dip into the integral 
equation , and apart from the diagonal term the only 
non-zero integral is given by 



Second term in (jClsp — The second 
term in the right hand side of Eq. (jClSP 
can be calculated by noting that, due to ro- 
tational invariance, {vo\dpAI[E{0),0]\vo) = 
lE,=x,y,zi'>'o\dl^M[EiO),0]\vo).' The kernel of the 
operator M'^^^[E{0),0] = ^M[£;(0), 0]/3 is then 



d^k' 



(2^)3 Fk'kq(i^, 0) ' 



(D2) 



where E — A_Edi,„(0) is the dimeronic energy at P = 0. 
In spherical coordinates {k',9',(j>') this integral can be 
easily solved by performing the substitution (jC22l) . and 
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by integrating over the azimuthal angle [63| we arrive at 
the following integral equation 



1 

4^ 



f{k',q,x') 



k A q 



= 0, (D3) 



where as usual x' — cos 9'. 

For an odd function of total angular momentum £ = 1 
we write the ansatz 

T]i^q= ^Uy{k, q, 9) + '^Vy{k,q, 9), (D4) 

where 9 is the angle between k and q. We have already 
encountered in the Appendix[C]the procedure of inserting 



this ansatz into the integral equation (ITSl) . which leads 
to an integral equation for the dimeronic energy E = 
AEdim(O) at P = given by 

'^{Uk + Vk)[k,q,9,E] 

+ j{U, + V,)[k,q,9,E]^0, (D5) 

where Uk, Vk, Uq and Vq are given in Eqs. 
(|C28IC29IC30IC31I) . We then get the integral equation 
in its final form by setting to zero the two lines of Eq. 
(jPSP separately since they are linearly independent |6J]: 



Uk[k,q,9,E] + Vk[k,q,9,E]^Q 
Uq[k,q,9,E] + Vq[k,q,9,E]=Q. 
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